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ABSTRACT 


Silicon  based  semiconductor  devices  are  rapidly  approaehing  the  theoretieal  limit 
of  operation  and  are  beeoming  unsuitable  for  future  military  requirements.  The  scope  of 
semieonduetor  deviees  has  been  expanded  by  wide  bandgap  deviees  such  as  gallium 
nitride  (GaN)  to  inelude  the  possibility  for  high  power  and  high  frequency  operation.  A 
new  generation  of  high  speed  -  high  frequeney  deviees  is  required  to  meet  eurrent  and 
future  military  needs.  The  Gallium  Nitride  High  Eleetron  Mobility  Transistor  (HEMT)  is 
showing  great  promise  as  the  enabling  teehnology  in  the  development  of  military  radar 
systems,  eleetronie  surveillanee  systems,  eommunieations  systems  and  high  voltage 
power  systems.  Typieally,  sapphire  or  silieon  earbide  is  utilized  as  the  substrate  material 
in  most  HEMT  designs.  This  thesis  explores  the  possibility  of  utilizing  a  diamond 
substrate  to  inerease  the  power  handling  eapability  of  the  AlGaN/GaN  HEMT.  Diamond 
offers  inereased  thermal  property  parameters  that  can  be  simulated  in  the  commereially 
available  Silvaeo  software  paekage.  A  eomplete  eleetrieal  and  thermal  analysis  of  the 
model  was  eondueted  and  eompared  to  aetual  deviee  oharaeteristies.  The  results  of  the 
software  simulation  and  measurements  on  the  test  deviees  indieate  diamond  substrates 
will  enable  the  HEMT  to  be  operated  at  a  higher  power  than  traditional  sapphire  substrate 
HEMTS. 
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EXECUTIVE  SUMMARY 


Current  AlGaN/GaN  High  Electron  Mobility  Transistors  (HEMTs)  grown  on 
sapphire  substrates  have  a  low  reliability  due  to  thermal  effects.  The  need  for  a  reliable 
high-frequency,  high-power  HEMT  device  is  necessary  for  future  naval  communication 
and  radar  needs.  A  computer  model  of  an  AlGaN/GaN  HEMT  has  been  developed  for 
the  purpose  of  studying  the  thermal  effects  due  to  self  heating.  The  thermal  effects  were 
examined  on  a  device  grown  on  a  sapphire  substrate  and  then  on  a  diamond  substrate  to 
determine  if  the  high  thermal  conductivity  parameter  of  diamond  would  provide  a 
reduced  device  operating  temperature.  The  model  developed  in  this  thesis  was  created 
with  the  commercially  available  Silvaco  software  package  and  includes  the  BLAZE™ 
heterojunction  routine,  the  GIGA^^  thermal  modeling  routine,  as  well  as  various  models 
to  simulate  the  low  field  mobility  as  a  function  of  doping  and  lattice  temperature. 

The  ATLAS'’’’^  model  was  created  based  on  the  parameters  extracted  from  an 
actual  AlGaN/GaN  HEMT  grown  on  a  sapphire  substrate.  The  parameters  extracted 
from  the  actual  device  were  utilized  as  baseline  for  the  model  development. 

One  of  the  critical  parameters  necessary  to  simulate  heterojunction  devices  is 
polarization.  The  polarization  effects  present  in  heterojunction  devices  are  spontaneous 
polarization  and  the  piezoelectric  effect  and  are  crucial  in  the  formation  of  the  unique  two 
dimensional  electron  gas  (2DEG)  channel  present  in  heterojunction  devices.  Although 
ATLAS™  provides  models  to  simulate  the  effects  individually,  the  results  were  less  than 
desirable.  In  the  NFS  model,  an  interface  charge  was  placed  on  the  AlGaN  side  of  the 
heterojunction  to  simulate  the  polarization  effects.  This  interface  charge  provided  the 
necessary  electron  concentration  to  form  the  2DEG  which  forms  an  abrupt  hangup  shift  at 
the  interface.  This  abrupt  change  forms  a  quantum  well  which  induces  the  2DEG  to 
form.  All  parameters  extracted  from  the  physical  model  were  able  to  be  re-created  in  the 
software  model  to  a  high  degree  of  accuracy. 

Simulations  were  conducted  utilizing  a  sapphire  substrate  to  match  the  operational 
parameters  of  the  physical  device.  Once  an  accurate  model  was  developed,  the  substrate 
material  was  changed  to  diamond,  and  the  thermal  characteristics  were  studied.  A 


XV 


simulation  of  the  sapphire  substrate  deviee  and  the  diamond  substrate  deviee  was 
eondueted  to  determine  if  a  ehange  in  the  thermal  parameters  affeeted  the  deviee  heating. 
The  sapphire  substrate  HEMT  model  was  validated  with  measurements  over  a 
temperature  range  of  21^C  (room  temperature)  to  300°C.  The  validated  deviee  modeling 
then  showed  a  dramatie  ehange  in  the  lattiee  temperature  as  well  as  the  ehannel 
temperature  of  the  deviee  when  the  substrate  was  ehanged  to  diamond.  The  diamond 
substrate  redueed  the  peak  temperature  of  the  deviee  by  1/2  while  the  average  ehannel 
temperature  was  redueed  by  1/3. 

The  work  presented  in  this  thesis  represents  the  most  aeeurate  NFS  HEMT  model 
developed  to  date  and  ean  be  utilized  in  follow-on  efforts  to  study  reliability  and  high 
frequeney  operations.  AlGaN/GaN  HEMT  is  showing  great  promise  in  beeoming  the 
enabling  teehnology  in  high-power,  high-frequeney  semieonduetor  amplifiers,  and  the 
results  of  this  thesis  indieates  that  utilizing  a  diamond  substrate  will  enhanee  the 
operation  of  the  HEMT. 
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I.  INTRODUCTION 


The  purpose  of  this  thesis  was  to  ereate  an  aeeurate  eleetrieal  and  thermal 
eomputer  model  of  an  AlGaN/GaN  HEMT  deviee  in  order  to  investigate  deviee  operating 
parameters  as  the  substrate  thiekness  and  material  are  ehanged.  To  develop  an  accurate 
model,  an  actual  AlGaN/GaN  HEMT  device  was  examined  and  parameters  extracted 
utilizing  Silvaco’s  EITMOST™  program.  The  biggest  departure  of  this  thesis  compared 
to  previous  work  in  this  area  is  that  this  work  was  based  and  verified  with  physical 
devices  and  correlated  with  simulated  models.  The  data  extracted  from  the  physical 
device  served  as  a  benchmark  in  the  development  of  the  device  model.  The  AlGaN/GaN 
HEMT  model  was  created  utilizing  Silvaco’s  ATEAS™,  BLAZE™^  and  GIGA™ 
routines. 

This  chapter  will  discuss  the  background  that  motivated  this  thesis,  current  uses  of 
HEMT  devices,  the  problems  associated  with  HEMT  devices,  ongoing  research  and 
funding  associated  with  HEMT  research,  and  will  conclude  with  a  summary  of  previous 
results  and  problems  encountered.  Chapter  two  will  describe  HEMT  device  physics  and 
operation.  Chapter  three  will  outline  the  procedures  for  testing  the  actual  HEMT  device, 
the  Silvaco  model  design  as  well  as  the  physics  models  utilized  to  match  the  model  with 
the  actual  device,  and  compare  the  physical  device  characteristics  with  the  model  results. 
Einally,  chapter  four  will  conclude  the  thesis  with  conclusions  and  recommendations  and 
discuss  possible  future  research. 

A,  BACKGROUND 

To  achieve  the  vision  of  Sea  Power  21,  the  United  States  Navy  requires  new 
technology  in  the  field  of  high-power,  high-frequency  electronics  to  field  new  state  of  the 
art  systems.  Current  silicon  technology  is  not  capable  of  achieving  the  performance 
required  by  future  systems.  The  technology  that  is  poised  to  fulfill  the  future  needs  of 
industry  and  the  Navy  is  wide  bandgap  (WBG)  semiconductor  devices.  In  particular,  the 
most  promising  device  to  satisfy  many  of  the  future  needs  of  the  Navy  is  the  GaN  High 
Electron  Mobility  Transistor  (HEMT).  Some  desirable  attributes  of  next  generation 
WBG  electronics  include  the  ability  to  withstand  currents  in  excess  of  5  kA  and  voltages 
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in  excess  of  50  kV,  provide  rapid  switching,  maintain  good  thermal  stability  while 
operating  at  temperatures  above  250°  C,  have  small  size  and  light  weight,  and  be  able  to 
function  without  bulky  heat-dissipating  systems  [1].  Table  1  provides  a  comparison  of 
the  various  semiconductor  technologies  available  today. 


Si 

(--) 

GaAs 
(AIGaAs/ 
In  GaAs) 

InP 

(InAIAs/ 

InGaAs) 

4H 

SIC 

(--) 

GaN 

(AIGaN/ 

GaN) 

Bandqap (eV) 

1.1 

1.42 

1.35 

3.26 

3.49 

Electron  mobility 
(cm^/Vs) 

1500 

8500 

10000 

700 

900 

Saturated  (peak) 
electron  velocity 
(x10^  cm/s) 

1 

2.1 

2.3 

2 

2.7 

2DEG  sheet  electron 
density  (cm'^) 

NA 

<4x10'^ 

<4x10'^ 

NA 

20x10’^ 

Critical  breakdown 
field  (MV/cm) 

0.3 

0.4 

0.5 

2 

3.3 

Thermal  conductivity 
(W/cm-K) 

1.5 

0.5 

0.7 

4.5 

>1.7 

Relative  dielectric 
constant  (fr) 

11.8 

12.8 

12.5 

10 

9.0 

Table  1:  Electronic  Properties  of  Semiconductor  Materials  [2]. 


In  addition  to  becoming  the  enabling  technology  for  future  systems,  GaN  HEMTs 
have  the  potential  to  replace  vacuum  tubes  in  high  power  millimeter  wave  radar  and 
electronic  warfare  systems.  They  have  the  potential  to  deliver  high  power  at  high  voltage 
and  do  it  at  a  fraction  of  the  cost  of  current  vacuum  tube  technology.  The  wide  bandgap 
of  GaN  increases  the  breakdown  field  by  live  times  and  the  power  density  by  a  factor  of 
10  to  20  compared  with  GaAs  based  devices.  This  results  in  smaller  components  with 
less  capacitance  for  the  same  operating  power,  which  results  in  a  wider  range  of 
bandwidth. 

GaN  HEMT  devices  appear  to  be  a  great  technology  on  paper,  and  great  progress 
is  being  made  toward  technology  development  and  process  improvement.  There  have 
been  demonstrations  of  very  high  output  power  densities  in  GaN  HEMTS,  but  significant 
developmental  work  remains.  One  of  the  remaining  issues  is  reliability.  The  output 
power  from  these  devices  can  permanently  degrade  over  a  relatively  short  period  of  time. 
A  phenomenon  known  as  current  collapse,  a  reduction  in  the  dc  drain  current,  after  the 

application  of  a  high  drain-source  voltage  is  prevalent  in  many  GaN  HEMT  devices. 
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The  high  voltage  leads  to  the  injeetion  of  hot  carriers  into  regions  of  the  device  adjacent 
to  the  conduction  channel  that  contains  deep  traps  [3].  One  reliability  issue  this  thesis 
attempts  to  address  is  how  increased  thermal  stress  affects  the  operation  of  the  GaN 
HEMT.  It  has  been  shown  that  the  DC  characteristics  are  better  for  HEMTs  fabricated  on 
SiC  substrate  than  on  a  sapphire  substrate,  possibly  because  of  the  increased  thermal 
conductivity  of  SiC.  The  model  developed  in  this  thesis  will  examine  the  DC 
characteristics  of  the  HEMT  utilizing  a  diamond  substrate.  Diamond’s  thermal 
conductivity  is  about  twice  that  of  SiC  and,  therefore,  will  reduce  the  channel 
temperature  of  the  HEMT,  which  will  increase  the  power  density  potential  and  reduce 
thermally  induced  reliability  issues. 

B,  ALGAN/GAN  HEMT  RESEARCH 

1.  Defense  Advanced  Research  Projects  Agency  (DARPA)/Microsystems 
Technology  Office  (MTO) 

DARPA  is  currently  funding  many  research  projects  dealing  with  WBG 
development  and  implementation.  One  DARPA  program  currently  underway  is  the  Wide 
Bandgap  Semiconductor-Radio  Erequency  (WBGS-RE)  program.  The  program 
objectives  are: 

•  Develop  reliable,  high  performance  devices  and  MMICs  with 

o  Higher  power 
o  Superior  thermal  performance 
o  Higher  efficiency 

•  Rapid  insertion  into  DoD  RE  systems. 

The  WBGS-RE  program  is  organized  into  three  phases.  Phase  I  is  complete  and  was 
focused  on  WBG  semiconductor  materials.  The  output  of  Phase  I  results  in  high  quality 
Si  and  SiC  substrates  and  highly  uniform  GaN  epitaxy,  utilizing  both  MBE  and  MOCVD. 
Phase  II  is  scheduled  for  three  years  and  is  currently  underway.  Phase  II  is  concerned 
with  device  technologies  such  as  X  and  Q  band  discrete  components.  The  output  of 
Phase  II  will  be  device  demonstrations.  Phase  III  is  a  two-year  phase  concerned  with 
electronic  integration  and  circuit  technologies. 
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The  DARPA  WBGS-RF  program  has  been  awarding  contracts  in  the  fields  of 
wide  bandgap  materials  optimization,  wide  bandgap  device  fabrication  and  modeling, 
high  yield  wide  bandgap  MMICs,  and  packaging  and  thermal  management  of  wide 
bandgap  devices.  There  were  multiple  contract  awards  starting  in  calendar  year  2005 
totaling  $75-$l  10  million  dollars. 

2.  DARPA  Wide  Bandgap  Semiconductor  Technology  Initiative 
(WBGSTI) 

Another  project  sponsored  by  DARPA  that  is  dedicated  to  advancing  WBG 
technology  is  the  WBGSTI.  The  objectives  of  the  WBGSTI  program  are  to  scale  up  high 
quality  SiC  substrates,  develop  alternative  substrates,  develop  uniform  AlGaN/GaN 
HEMTs,  and  to  examine  materials/device  correlations.  This  project  is  unique  because  it 
has  developed  a  Tri-Service  team  comprised  of  the  Air  Force  Research  Laboratory 
(AFRL),  Army  Research  Laboratory  (ARE),  and  the  Naval  Research  Laboratory  (NRL). 
The  Tri-Service  team  is  composed  of  a  Materials  Panel  and  a  Device  Panel.  Figure  1 
depicts  the  interaction  of  the  Tri-Service  team. 


4 


Figure  1:  Tri-Service  Team  Interaction  [4], 


The  WBGSTI  is  a  multi-phase,  multi-year  project  that  is  poised  to  revolutionize 
the  development  of  WBG  devices.  This  project  has  the  potential  to  change  the  way  in 
which  high-power,  high-frequency  systems  are  designed. 


C.  SUMMARY  OF  PREVIOUS  RESEARCH 

This  section  summarizes  previous  Naval  Postgraduate  School  modeling  attempts 
relating  to  AlGaN/GaN  HEMT  modeling  research.  The  previous  efforts  by  Eimers  [5], 
Holmes  [6],  and  Salm  [7]  concentrated  on  correctly  modeling  the  HEMT  device. 
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Previous  versions  of  the  Silvaeo  software  did  not  provide  adequate  models  for  the  surfaee 
eharge  or  piezoeleetrie  effeets  eontained  in  AlGaN/GaN  HEMTS.  Therefore,  many  of 
the  earlier  attempts  to  eorreetly  model  the  HEMT  deviee  were  foeused  on  independently 
developing  these  models  to  inelude  in  the  Silvaeo  routine. 

In  [6],  Holmes  attempted  to  ineorporate  the  polarization,  piezoeleetrie  and 
spontaneous  effeets  into  the  HEMT  model  utilizing  a  C-interpreter  method.  Holmes  was 
able  to  get  an  adequate  I-V  eharaeteristie  of  a  HEMT  deviee  utilizing  a  C-interpreter 
routine.  However,  his  model  needed  to  be  refined  to  more  elosely  mateh  the  overall 
operation  of  the  HEMT  deviee. 

Sinee  the  eompletion  of  Holmes’  thesis  in  2002,  Jon  C.  Ereeman  [8]  developed  a 
method  to  model  GaN  HEMTs  with  the  Silvaeo  software  paekage.  In  partieular,  he 
utilized  ATEAS™  to  add  the  effeets  of  strain,  piezoeleetrieity,  and  spontaneous 
polarization  at  the  interfaees.  The  objeetives  of  Holmes’  thesis  seem  to  have  been  met  by 
Silvaeo’s  software  upgrade  and  Ereeman’s  researeh.  The  model  developed  in  this  thesis 
ineorporated  the  new  Silvaeo  models  and  resulted  in  a  more  aeeurate  deviee  model. 

Salm’s  [7]  thesis,  was  the  first  Naval  Postgraduate  Sehool  attempt  to  utilize  the 
piezoeleetrie  model  developed  by  Silvaeo  in  a  HEMT  deviee  simulation.  Although  he 
was  not  able  to  get  aeeurate  results  with  the  piezoeleetrie  model,  he  modeled  the  effeet 
with  an  interfaee  eharge,  and  his  results  were  adequate  and  more  aeeurately  modeled  the 
desired  eharaeteristies  of  a  HEMT  deviee.  His  thesis  also  studied  the  thermal 
eharaeteristies  of  the  deviee  utilizing  the  GIGA^m  module.  The  GIGA™  module  is 
eapable  of  giving  a  numerieal,  as  well  as  illustrative,  representation  of  deviee  heating. 
Salm’s  results  indieated  a  remarkable  differenee  in  deviee  heating  due  to  a  ehange  in 
substrate  material.  The  speeifies  of  the  strueture  and  deviee  physies  will  be  diseussed  in 
ehapter  three  of  this  thesis. 
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II.  HIGH  ELECTRON  MOBILITY  TRANSISTORS  (HEMT) 


This  chapter  presents  an  introduction  to  the  High  Eleetron  Mobility  Transistor 
(HEMT).  The  theory  of  operation  for  a  HEMT  device  will  be  discussed  to  inelude 
equations,  bandgap  ealculations,  and  the  unique  feature  of  the  HEMT,  the  two 
dimensional  eleetron  gas  (2-DEG).  This  chapter  will  conclude  with  an  overview  of  the 
AlGaN/GaN  structure  designed  and  fabricated  by  the  Elniversity  of  California  at 
Berkeley. 


A,  HEMT  FUNDAMENTALS 

1,  HEMT  Structure 

The  GaN  HEMT  is  similar  to  a  typieal  MESEET  in  that  it  is  a  three  terminal 
device  that  utilizes  the  gate  to  control  current  flow.  The  Naval  Postgraduate  School  was 
provided  with  an  array  of  devices  designed  and  tested  by  Susie  Tzeng  [9]  from  the 
University  of  California  at  Berkeley.  The  exact  structure  dimensions  shown  in  Eigure  3 
selected  for  this  thesis  includes  a  gate  length  (Eg)  of  Ipm,  a  gate  width  (Wg)  of  100  pm, 
and  a  source  to  gate  (Egs)  and  gate  to  drain  (Egd)  distance  of  2  pm.  Additionally,  the  A1 
(x)  composition  of  the  Berkeley  HEMT  was  varied  between  28%  and  35%.  The  A1 
composition  selected  for  this  thesis  was  28%.  Eigure  2  illustrates  the  layer  structure  of 
the  device  and  Eigure  3  defines  the  regions  for  the  dimensions  described  above. 


1.5um  GaN 


Figure  2:  AlGaN/GaN  HEMT  Layer  Structure, 
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Figure  3:  HEMT  Device  Dimensions  [10], 

Ti/Al/Ni/Au  ohmic  metallization  was  deposited  by  e-beam  evaporation  and 
patterned  by  lift-off  [9].  This  proeess  resulted  in  a  eontaet  resistance  of 
approximately  5  X  lO^^Q-cm^ .  The  physical  characteristics  of  the  Berkeley  device  were 
utilized  as  an  input  to  the  ATLAS'^’^  model  and  are  diseussed  in  chapter  three  of  this 
thesis.  Figure  4  is  a  scanning  electron  microscope  (SEM)  image  of  the  HEMT  deviee 
structure  where  the  gate  eontaet  ean  be  seen  positioned  between  the  source  and  drain 
ohmic  contacts  [9]. 

Although  a  complete  strain  analysis  of  the  deviee  is  beyond  the  seope  of  this 
thesis,  an  overview  of  the  strains  and  polarization  effeets  in  a  HEMT  deviee  will  be 
diseussed  in  seetion  two  of  this  ehapter.  The  polarization  effeets  inherent  in  HEMT 
devices  have  been  modeled  by  Dr.  Jon  Ereeman  [8]  and  ean  be  ineorporated  utilizing  the 
Silvaeo  software  paekage.  The  modeling  efforts  of  this  thesis  will  be  detailed  in  ehapter 
three  of  this  thesis  to  inelude  the  polarization  effeets  that  occur  in  HEMT  devices. 
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Figure  4:  SEM  Image  of  Berkeley  GAN  HEMT  device.  [9], 

2,  AlGaN/GaN  HEMT  Fundamentals 

The  AlGaN/GaN  HEMT  also  known  as  the  modulation  doped  field  effect 
transistor  (MODFET)  has  similar  DC  characteristics  to  the  MESFET.  The  unique  feature 
of  the  HEMT  is  channel  formation  from  carriers  accumulated  along  a  grossly  asymmetric 
heterojunction,  i.e.  a  junction  between  a  heavily  doped  high  bandgap  and  a  lightly  doped 
low  bandgap  region  [1],  Instead  of  the  channel  being  formed  as  in  a  standard 
homojunction  device,  a  two  dimensional  electron  gas  (2DEG)  is  formed  at  the 
heterointerface.  To  achieve  proper  operation  of  the  device,  the  barrier  layer  AfiGai-xN 
must  be  at  a  higher  energy  level  than  the  conduction  band  of  the  GaN  channel  layer.  This 
conduction  band  offset  transfers  electrons  from  the  barrier  layer  to  the  channel  layer.  The 
electrons  that  are  transferred  are  confined  to  a  small  region  in  the  channel  layer  near  the 
heterointerface.  This  layer  is  called  the  2DEG  and  a  defining  characteristic  of  the 
HEMT.  There  are  many  factors  that  determine  the  quality  of  the  2DEG.  The  factors 
involved  in  the  development  of  the  2DEG  are  type  of  substrate,  growing  method,  and 
level  of  doping  of  the  carrier  supply  layer  [10].  Electron  density  and  mobility  are 
important  factors  when  dealing  with  a  2DEG  and  a  more  detailed  discussion  of  the  2DEG 
will  be  presented  later  in  this  chapter. 
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a.  HEMT  Equations 

This  section  will  outline  basic  HEMT  parameters.  The  primary  operating 
parameters  dealing  with  high  power  HEMT  operations  are  the  transconductance  (gm),  the 
2DEG  sheet  charge  concentration  (ns),  the  drain  current  (Id),  the  threshold  voltage  (Vt), 
and  the  contact  barrier  height  ( ^4,  )• 

The  typical  DC  measurement  to  determine  the  operation  of  a  HEMT  is  the 
drain  current  versus  the  voltage  drain  to  source  (Id  vs.  Vds)  while  the  gate  is  swept  from 
below  cutoff  to  saturation.  Eigure  5  is  an  Id  vs.  Vds  curve  extracted  from  one  of  the 
Berkeley  devices. 


Room  Temperature  IDS  Vs.  VDS  curve 
Data  from  5_12_0C.ssf 


Figure  5:  DC  Output  characteristic  of  AlGaN/GaN  HEMT, 
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Unlike  traditional  MOS  devices,  the  charge  flowing  in  a  HEMT  device  is 
normally  conflned  to  the  2DEG  and  changed  by  the  gate  voltage.  When  a  gate  voltage  is 
applied  to  the  device,  the  sheet  charge  concentration  (ns)  is  changed  proportionally  to 
gate  voltage  and  is  given  by  equation  2.1  [11]: 


n,  =  ■ 


q{d^  +  AJ) 


(2.1) 


where  and  are  the  dielectric  permeability  and  the  thickness  of  the  wide  bandgap 

semiconductor,  respectively,  and  AJ  can  be  interpreted  as  the  effective  thickness  of  the 
2DEG  [11].  A  typical  ns  value  for  an  AlGaN/GaN  heterojunction  is  ~  lO'^/cm^. 

The  threshold  voltage  is  an  important  parameter  because  it  is  a  measure  of 
when  the  HEMT  device  will  begin  conducting.  Vx  can  be  expressed  as: 

—  d"  (2.2) 

la, 


where  is  the  gate  Schottky-barrier  height,  and  AE^  is  the  change  in  the  conduction 
band  at  the  heterojunction,  and  Nd  is  the  background  doping  of  the  GaN  layer. 

Saturation  current  and  transconductance  are  important  measures  of  a 
HEMT  device  when  utilized  as  a  power  amplifier.  Eor  long  channel  devices  (Eg> 
0.25pm),  such  as  the  ones  presented  in  this  thesis,  losat  can  be  expressed  the  same  as  a 
standard  MOSEET  [9]: 


/ 


DSat 


!^CW 

IL 


(K-Kf 


(2.3) 


where  //  is  the  carrier  mobility,  C  is  the  gate  capacitance,  W  is  the  gate  width,  and  E  is 
the  gate  length. 

One  of  the  measures  of  quality  for  a  HEMT  device  is  transconductance.  It 
is  a  measure  of  how  the  drain  current  will  change  with  a  change  in  the  gate  voltage. 
Mathematically,  transconductance  is  defined  as  [12]: 
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The  output  conductance  plays  a  crucial  role  in  determining  the  output  impedance 
matching  properties  of  different  stages  of  an  amplifier. 

The  equations  presented  above  are  a  good  starting  point  for  modeling  the 
Berkeley  HEMT  device.  However,  there  are  many  variables  within  the  equations  that 
will  be  modified  to  achieve  a  proper  model  of  the  device.  There  are  a  variety  of  models 
included  in  the  Silvaco  package  that  can  be  incorporated  into  the  device  definition,  and 
when  necessary  a  more  detailed  description  of  the  chosen  model  will  be  provided. 

b.  HEMT  Material  Composition  Issues 

As  mentioned  above,  the  HEMT  is  a  three  terminal  device  where 
conduction  is  controlled  by  the  gate  voltage.  The  biggest  operational  difference  between 
a  HEMT  and  a  traditional  MOS  device  is  that  instead  of  creating  a  channel  between  the 
source  and  drain  with  a  doped  material,  a  wider  bandgap  material  barrier  layer 
(AlxGai-xN)  is  grown  on  top  of  the  channel  and  serves  as  a  carrier  supply  layer  [9]. 
Therefore,  the  channel  is  formed  as  a  2DEG  below  the  doped  layer  which  increases 
mobility  due  to  a  reduction  in  impurity  scattering.  A  technique  utilized  to  increase  the 
electron  density  in  the  channel  is  to  increase  the  A1  mole  fraction  (x)  in  the  AfiGai-xN 
layer.  An  increase  in  the  A1  mole  fraction  results  in  a  higher  bandgap  and  a  larger 
conduction  band  discontinuity  AE^  relative  to  the  GaN  [13].  The  higher  bandgap  and 

AE^,  results  in  an  increase  in  the  composite  breakdown  field  and  improves  carrier 

confinement.  Additionally,  the  piezoelectric  and  spontaneous  polarization  charges 
increase  and  serve  to  confine  the  carriers  to  the  2DEG.  The  effects  of  increasing  the  A1 
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mole  concentration  are  attractive,  but  there  is  a  limit  to  benefit  gained.  Figure  6 
illustrates  the  electron  density  as  a  function  of  A1  composition. 


Figure  6:  2DEG  electron  density  distribution  for  different  Al(x)  composition 

[14]. 

Increasing  the  A1  mole  fraction  also  results  in  a  reduced  carrier  mobility  due  to  different 
effects  including,  intersubband  scattering,  increased  alloy  disordering  in  AlGaN, 
increased  density  of  interface  charges  and  larger  potential  fluctuations  at  the  interface  due 
to  the  surface  roughness  [14].  The  A1  composition  for  the  device  tested  in  this  thesis  is 
28%,  which  roughly  correlates  to  an  electron  density  of  2el9  cm'  . 
c.  HEMT  Operational  Parameters 

This  section  will  provide  an  overview  of  the  device  structure,  where  the 
2DEG  is  created,  and  the  energy  bandgap  diagrams  of  an  AlGaN/GaN  HEMT  device. 

The  HEMT  operational  characteristics  are  based  on  the  MESEET.  The 
2DEG  is  formed  just  below  the  heterojunction  where  the  defects  are  lower  and  a  higher 
electron  mobility  is  possible.  Eigure  3  gives  a  visual  estimate  of  where  the  2DEG  will  be 
within  the  structure.  The  conductivity  of  the  two  dimensional  channel  is  given  by  [10]; 

cr  =  qn^/u  (2.6) 

where  q  is  the  electron  charge,  is  the  sheet  carrier  concentration  given  by  equation  2.1, 
and  jj.  is  electron  mobility. 
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Figure  7  depicts  the  typical  bandgap  diagram  of  an  AlGaN/GaN  HEMT 

device. 


Figure  7:  Bandgap  diagram  of  AlGaN/GaN  HEMT  [9], 

One  of  the  defining  characteristic  of  a  heterojunction  device  is  the  abrupt  bandgap  shift 
where  the  AlGaN  layer  meets  the  GaN  Layer.  The  polarization  effects  that  occur  at  this 
interface  allows  the  2DEG  to  be  formed  and  become  a  conducting  channel  when  properly 
biased. 

The  Id  vs.  Vds  curve  generated  in  Eigure  5  will  serve  as  the  basis  for 
discussion  of  the  HEMT  operating  conditions.  The  initial  values  for  the  measurement 
were: 

•  Source  voltage  (Vs)  =0 

•  Drain  voltage  (Vd)  =  0 

•  Gate  voltage(Vg)=-5 

The  source  voltage  was  held  constant  and  the  drain  voltage  was  swept  from  0-7  V.  The 
gate  voltage  began  at  -5  V  and  was  incremented  by  1  V  per  sweep  to  a  final  value  of -IV. 
With  the  Vs  and  Vg  =  OV,  Vd  was  incremented  in  steps  up  to  7  V.  It  is  apparent  from 
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Figure  5  that  a  small  current  began  to  flow  in  the  2DEG  channel  with  no  gate  voltage. 
By  applying  a  positive  voltage  to  the  drain,  current  will  flow  due  to  the  potential 
difference  between  the  source  and  drain,  but  the  magnitude  of  the  current  flow  is 
controlled  by  the  gate.  Each  branch  in  the  diagram  represents  a  1  V  increase  in  Vg  which 
causes  the  drain-source  current  to  increase  linearly  until  the  saturation  point  of  the  device 
is  reached.  The  value  of  losat  is  dependant  upon  the  concentration.  By  increasing  the 

concentration  of  ,  the  2DEG  is  capable  of  transferring  more  charge  from  source  to 
drain. 

Unlike  a  typical  MOS  U  vs.  Vds  curve,  GaN  HEMTs  on  a  sapphire 
substrate  show  a  small  negative  slope  in  the  saturation  region.  This  anomaly  is  known  as 
negative  resistance  and  is  attributed  to  self-heating  of  the  device  due  to  the  poor  thermal 
conductivity  of  sapphire[l]. 

3.  Two  Dimensional  Electron  Gas  (2DEG) 

The  2DEG  in  AlxGai.xN/GaN  based  HEMT  is  induced  by  strong  polarization 
effects  [15].  Ambacher  et  al.,  utilizing  piezoelectric  constants  obtained  from  previous 
research,  calculated  the  induced  sheet  charge  bound  at  the  heterojunction  interface.  His 
conclusion  was  that  the  formation  of  the  2DEG  in  AlGaN/GaN  structures  was  the  result 
of  both  spontaneous  polarization  (Psp)  and  the  strain  induced  piezoelectric  effects  (Ppe) 
[15].  The  Silvaco  software  package  includes  new  modeling  parameters  that  will  enable  a 
more  realistic  HEMT  device  model  to  be  developed.  The  following  sub-sections  address 
the  polarization  effects  of  a  GaN  HEMT  and  will  provide  some  of  the  variables  that 
affect  the  2DEG. 

a.  Spontaneous  Polarization 

The  two  distinct  polarization  mechanisms  found  in  GaN  HEMTs  enable  a 
high  carrier  concentration  even  though  there  may  be  no  intentional  doping.  Spontaneous 
polarization  is  the  effect  that  is  present  at  the  heterojunction  even  though  there  is  no 
strain  present.  The  spontaneous  polarization  (Psp)  is  developed  due  to  the  wurtzite  group 
Ill-nitrides  being  tetrahedrally  coordinated  with  a  lack  of  symmetry  along  the  [0001]  axis 
or  c-direction  [10].  The  polarization  effects  vary  according  to  the  how  the  structure  was 
formed.  The  structure  will  have  either  N-face  or  Ga-face  surfaces.  The  Berkeley  device 
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provided  was  grown  on  sapphire  by  MOCVD  whieh  results  in  Ga-faeed  struetures  only. 
Figure  8  shows  the  erystal  strueture  of  a  wurtzite  Ga-faee. 


Gs-fsce 


Substrate 


Figure  8:  Schematic  drawing  of  the  crystal  structure  of  wurtzite  Ga-face[15], 

In  the  ease  of  the  wurtzite  erystal  grown  in  the  manner  previously 
diseussed,  the  Psp  will  be  given  as  a  eonstant  determined  by  the  method  of  growth.  Table 
3  provides  the  spontaneous  polarization,  piezoeleetrie  and  dieleetrie  eonstants  of  AIN  and 
GaN.  Ambaeher  et  al.  [15]  developed  this  table  from  many  sources  and  utilized  the 
constants  to  determine  the  Psp  and  PpE  for  an  AlGaN/GaN  structure. 


woutzite 

.41N 

GaN 

Psp  [C  m-] 

-0.081 

-0.029 

1.46 

0.73 

1.55 

1 

[Cm^] 

0.65 

0.44 

«3l 

-0.60 

-0.49 

[Cm-] 

-0.58 

-0.36 

0.33 

—0.22 

^15 

-0.48 

-0.3 

[Cmh 

-0.33 

—0.22 

«11 

9.0 

9.5 

^33 

10.7 

10.4 

Tahle  2:  Spontaneous  polarization,  piezoelectric  and  dielectric  constants  of 

AIN,  GaN  [15]. 

The  ability  to  achieve  a  high  carrier  concentration  at  the  2DEG  is  determined  by  the 
ability  of  the  structure  to  confine  carriers  through  the  polarization  effects.  To  achieve  a 
high  carrier  confinement,  the  polarization  directions  should  be  the  same.  The 
polarization  induced  sheet  charge  density  (cr)  should  be  positive.  This  positive  sheet 
charge  attracts  electrons  and  forms  the  2DEG.  Eigure  9  illustrates  the  polarization  charge 
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effects  for  the  different  types  of  strain.  The  first  case  (a)  where  AlGaN  (relaxed)  is 
grown  on  top  of  GaN  (relaxed),  a  positive  sheet  charge  ( +cr )  is  developed  at  the  interface 
and  a  2DEG  forms.  To  enhance  the  polarization  effects,  example  (b)  grew  tensile 
strained  AlGaN  over  GaN  which  resulted  in  all  polarization  effects  acting  in  the  same 
direction.  This  also  results  in  a  positive  sheet  charge  at  the  interface.  The  piezoelectric 
(Ppe)  effects  will  be  illustrated  and  calculated  in  the  next  section. 


a) 


Ga-face 


AlGaN 

relaxed 


GaN 

relaxed 


AlCeN 

tensile  strain 


GaN 

relaxed 


GaN 

compressive 

strain 


AlGaN 

relaxed 


Figure  9:  Ga-Faced  polarization  induced  sheet  charge  density  and  directions  for 
the  spontaneous  and  piezoelectric  polarization  effects  [15]. 


b.  Piezoelectric  Polarization 

The  total  polarization  in  the  structure  is  given  by  Psp  +  Ppe=Ptotal- 
Utilizing  Table  2,  the  piezoelectric  polarization  can  be  calculated  by  the  following 
equation  given  by  [15]; 

^PE  ~  Gl  (g  )  C^-^) 
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where  ao  and  Co  are  the  equilibrium  values  of  the  lattice  parameters, 


V  <^0  y 


is  the 


strain  along  the  c-axis,  and  the  in-plane  strain  a^=  a  = 


^  a- 
V  y 


is  assumed  to  be 


isotropic.  The  third  independent  component  of  the  piezoelectric  tensor  is  related  to 


the  polarization  induced  by  shear  strain,  which  is  not  applicable  in  this  device.  The 
relation  between  the  lattice  constants  of  the  hexagonal  GaN  is  given  by  Ambacher  et  al. 
[15]  as: 


(c  \ 

*-^13 

0 

1 

0 

=  -2 

0 

1 

1  ^0  y 

f  Q3  y 

K  a,  ) 

(2.8) 


where  Cb  and  C33  are  elastic  constants  provided  in  previous  literature.  The  total 
piezoelectric  polarization  in  the  direction  of  the  c-axis  can  be  determined  by  equations 
2.7  and  2.8  [15]: 


P  =2 

^  PE  ^ 


y  y 

Q-Qq 

y 

'^13 

V  y 

Gi  “G3 
V 

v^33  jy 

(2.9) 


Since  the  value  of 


i^T 

Gi  G3 

1^33 

<  0  for  AlGaN  over  the  entire  range  of  compositions. 


the  piezoelectric  polarization  is  negative  for  tensile  and  positive  for  compressive  strained 
barriers,  respectively  [15].  GaN  based  semiconductors  are  attractive  in  HEMT  designs 
due  to  the  high  piezoelectric  constants.  The  piezoelectric  fields  developed  in  GaN 
devices  can  create  a  vertical  electric  field  in  the  MV/cm  range  which  increases  the 
channel  conductivity.  Many  of  the  polarization  effects  outlined  above  have  been  modeled 
into  ATLAS™  as  an  upgrade  to  heterojunction  modeling. 


B,  BERKELEY  HEMT  DESIGN 

All  HEMT  structures  provided  to  NFS  were  metal-organic  chemical  vapor 
deposition  (MOCVD)  grown  by  Emcore  Corporation  [9].  The  structure  was  grown  on 
[0001]  c-plane  sapphire  substrate. 
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Figure  2  depicts  the  layer  by  layer  structure  of  the  device  tested.  From  bottom  to 
top,  it  consists  of  the  sapphire  layer,  a  buffer  layer,  1.5  pm  of  unintentionally  doped 
(UID)  GaN,  and  approximately  25  nm  of  UID  AlxGai-xN.  The  A1  composition,  (x),  in  the 
Berkeley  devices  was  varied  between  28%  and  35%  to  assess  the  impact  of  A1 
composition  on  the  2DEG  transport  properties  and  the  device  performance  [9].  The  GaN 
HEMT  was  fabricated  by  a  3-step  process:  dry-etch  mesa  isolation,  ohmic  contact 
metallization,  and  gate  metallization.  Eigure  4  provides  a  SEM  image  of  the  device. 
Table  3  [9]  shows  the  AlxGai.xN  thickness  as  a  function  of  the  percent  of  Al. 


Sample 

Number 

Al% 

AlxGai.xN 

thickness  (A) 

1 

28 

267 

2 

29 

230 

3 

30.2 

230 

4 

32 

221 

5 

34 

242 

6 

35 

206 

Table  3  Al  Composition  and  layer  thickness  of  AEGai-xN  barrier  layer  in 

AlxGai  xN/GaN  HEMT  structure  [9]. 
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III.  AlGaN/GaN  HEMT  SILVACO  MODEL 


The  model  presented  in  this  thesis  was  generated  with  Silvaco  International’s 
software  package.  An  input  deck  was  generated  in  the  DECKBUILD™  VWF  interactive 
tool,  solved  through  the  ATLAS  routine,  and  analyzed  utilizing  the  TONYPLOT  tool. 
Additionally,  the  BLAZE™  routine,  a  routine  specifically  designed  for  group  III-V 
materials  and  devices  with  position  dependent  band  structures  was  utilized  to  simulate  the 
AlGaN/GaN  HEMT. 

The  desired  result  of  this  thesis  was  to  generate  an  accurate  HEMT  model  and 
determine  the  thermal  characteristics  of  the  device  as  the  substrate  thickness  and  material 

TM 

are  changed.  The  GIGA  routine  was  utilized  to  determine  general  thermal 
characteristics  of  the  device.  GIGA  extends  the  ATLAS  program  to  account  for 
lattice  heat  flow  and  general  thermal  environments.  The  routing  implements  Wachutka’s 
thermodynamically  model  of  lattice  heating  which  accounts  for  Joule  heating,  heating 
and  cooling  due  to  carrier  generation  and  recombination,  and  the  Peltier  and  Thomson 
effects  [16]. 

The  Silvaco  software  suite  is  a  powerful  simulation  software  that  uses  physics 
based  simulations.  Physically-based  device  simulators  predict  the  electrical 
characteristics  that  are  associated  with  specified  physical  structures  and  bias  conditions. 
ATLAS™  achieves  this  simulation  by  applying  a  two  or  three  dimensional  grid  that  forms 
nodes  throughout  the  structure.  The  nodes  are  solved  in  a  matrix  with  charge  transport 
equations  from  Maxwell’s  continuity  and  drift-diffusion  equations  [16]. 

The  ATLAS™  simulator  is  a  command  driven  simulator  that  requires  an  input  in  a 
specific  sequence.  The  order  in  which  the  statements  occur  is  important  to  the  proper 
simulation  of  a  device.  There  are  five  groups  of  statements,  and  if  not  correctly  stated 
will  result  in  errors  or  erroneous  simulations.  The  required  groups  and  statements  are  as 
follows  [16]; 
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•  Structure  specification 

o  Mesh 
o  Region 
o  Electrode 
o  Doping 

•  Material  Models  Speeification 

o  Material 
o  Models 
o  Contaet  Interface 

•  Numerieal  Method  Selection 

o  Method 

•  Solution  Specifieation 

o  Log 
o  Solve 
o  Load 
o  Save 

•  Results  Analysis 

o  Extract 
o  TONYPLOT. 

The  above  eommands  represent  the  entire  structure  definition  and  solution  method.  It  is 
critical  to  have  a  firm  understanding  of  the  simulation  variables  prior  to  defining  a 
structure.  A  small  change  in  one  variable  or  solution  method  ean  have  undesirable 
consequenees  on  device  performanee.  Ligure  10  depicts  a  flow  chart  of  the  device 
simulation  utilizing  Silvaco  International’s  simulators  and  modules. 
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(Runtime  Environment) 


Figure  10:Silvaco  Simulation  Flowchart  [16], 

A,  SILVACO  SEMICONDUCTOR  MODELING  EQUATIONS 

The  models  developed  to  simulate  the  operation  of  semiconductor  devices  consist 
of  equations  derived  from  Maxwell’s  laws,  Poisson’s  equation,  the  continuity  equations 
and  the  drift-diffusion  transport  equations  [16].  The  following  sections  discuss  the 
equations  in  broad  terms,  and  when  necessary,  a  more  detailed  explanation  of  the 
equations  will  be  presented. 

The  method  in  which  the  above  equations  will  be  solved  is  a  prime  consideration 
when  developing  a  device  model  in  ATLAS™.  Several  different  numerical  methods  can 
be  used  for  calculating  the  solutions  to  various  device  structures.  Different  solution 
methods  are  employed  depending  on  the  situation.  The  model  developed  in  this  thesis 
incorporated  the  lattice  heating  parameter  simulated  through  GIGA’^”.  When  the  lattice 
heating  model  is  added  to  drift-diffusion  equations  an  extra  equation  is  added.  The 
BLOCK  algorithm  solves  the  three  drift  diffusion  equations  as  a  Newton  solution  and 
follows  this  with  a  decoupled  solution  of  the  heat  flow  equations  [16].  The  NEWTON 
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algorithm  solves  all  four  equations  in  a  eoupled  manner.  NEWTON  is  preferred  onee  the 
temperature  is  high,  however  BLOCK  is  quieker  for  low  temperature  gradients.  The 
eombination  utilized  for  this  thesis  was:  METHOD  BLOCK  NEWTON. 

1.  Poisson’s  Equation 

Poisson’s  equation  is  a  well  known  partial  differential  equation  that  has  funetions 
in  eleetrostaties  and  theoretieal  physios.  In  semioonduotor  modeling  it  often  serves  as  the 
starting  point  in  obtaining  quantitative  solutions  for  eleotrostatio  variables  [12].  In 
eleetrostaties,  Poisson’s  equation  relates  the  eleotrostatic  potential  to  the  spaoe  oharge 
density  and  is  given  by  [16]: 

div(£V'F)  =  -/?  (3.1) 

where  'P  is  the  eleotrostatio  potential,  s  is  the  looal  permittivity,  and  p  is  the  looal  spaoe 
oharge  density.  The  looal  spaoe  charge  density  is  a  function  of  all  mobile  and  fixed 
charges  including  electrons,  holes,  and  impurities.  The  electric  field  is  obtained  from  the 
gradient  of  the  potential  and  is  given  by  [16]: 


^  =  V'P  (3.2) 

2.  Carrier  Continuity  Equations 

The  carrier  continuity  equations  for  electrons  and  holes  are  defined  by  [16]: 


dn 

dt 


UwJ  +G„-R, 

q 


(3.3) 


dp 

dt 


(3.4) 


hole  current  densities,  Gn  and  Gp  are  the  generation  rates  for  electrons  and  holes,  Rn  and 
Rp  are  the  recombination  rates  for  electrons  and  holes  and  q  is  the  magnitude  of  the 
charge  on  the  electron.  It  is  possible  to  solve  for  only  the  holes  or  electrons  in  the 
equations,  and  the  model  presented  in  this  thesis  solves  the  equations  for  electrons  only. 
Hole  transport  in  the  HEMT  is  insignificant. 
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3,  Transport  Equations 

Electrons  in  thermal  equilibrium  at  temperature  Tl  within  a  semiconductor  lattice 
obey  Fermi-Dirac  statistics.  That  is  the  probability  f(^)  that  an  available  electron  state 
with  energy  a  is  occupied  by  an  electron  is  given  by  [16]; 


1 


f 

1  +  exp 

V 


(3.5) 


where  Ep  is  a  spatially  independent  reference  energy  known  as  the  Fermi  level  and  k  is 
the  Boltzmann’s  constant.  If  ^  -  Ef»  kEp,  equation  4.5  can  be  approximated  as: 


r 

f{s)  =  exp 

V 


KTl  ) 


(3.6) 


Statistics  based  on  equation  3.6  are  know  as  the  Boltzmann  statistics  and  were 
utilized  in  the  formulation  of  the  HEMT  model.  The  current  density  equations  know  as 
the  charge  transport  models  are  obtained  by  applying  simplifications  to  the  Boltzmann 
transport  equation.  These  assumptions  can  result  in  a  number  of  different  transport 
models  such  as  the  drift-diffusion  model,  the  energy  balance  model  or  the  hydrodynamic 
model  [16].  The  choice  of  charge  transport  model  will  then  have  a  major  influence  on 
the  choice  of  generation  and  recombination  model.  The  drift  diffusion  transport  model 
was  utilized  in  the  development  of  this  model.  The  following  equations  were  derived  in 
chapter  three  of  [16],  and  show  the  conventional  form  of  the  drift-diffusion  equation; 


=qnM„E^+qDyn 


J p  =  qp/^pEp  +  qDpVp 


where  Dn  and  Dp,  in  the  case  of  the  Boltzmann  statistics  corresponds  to: 

n 

Dn 

q 


n 

q 


(3.7) 

(3.8) 


(3.9) 

(3.10) 
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4,  Polarization  Effects 

Perhaps  one  of  the  most  ehallenging  aspeets  of  modeling  heterojunetion  deviees 
in  ATLAS^^”  is  implementing  the  polarization  effeets  present  at  the  interfaee  of  the 
heterojunetion.  It  has  been  widely  published  [15,  17]  that  surfaee  donor  like  traps  are  the 
souree  of  eleetrons  in  the  ehannel  and  the  polarization  effeets  in  AlGaN/GaN  struetures 
foree  the  eleetrons  into  the  ehannel.  Previous  modeling  efforts  have  modeled  the 

f  e  ^ 

polarization  effeets  as  an  interfaee  eharge  of  approximate  lei 3  — j  at  the 

yam  J 

heterojunetion  interfaee. 

Reeent  upgrades  to  ATLAS™  inelude  the  ability  to  model  polarization  effeets 
direetly.  The  ATLAS™  polarization  model  was  designed  to  simulate  polarization  in 
wurtzite  materials  and  inelude  spontaneous  polarization  and  piezoeleetrie  polarization. 
The  model  allows  the  designer  to  input  the  piezoeleetrie  eonstants  esi  and  ess  as  well  as 
the  elastie  eonstants  Cis  and  Css.  Various  methods  to  implement  the  polarization  model 
into  the  deviee  model  were  attempted  however,  all  attempts  were  unsueeessful.  When 
ineorporated  into  the  model,  the  polarization  effeets  were  evident  on  the  eontour  plots, 
but  I  was  unable  to  get  a  representative  bandgap  diagram  or  Id  vs.  Vd  eharaeteristie. 

The  polarization  effeets  were  modeled  as  an  interfaee  eharge  of  approximately 

f  e  ) 

lel3  — Y  on  the  GaN  side  of  the  heterojunetion.  Ambaeher  et  al.  [15]  utilizes  the 
ycm  ) 

lathee  eonstant  (a),  the  piezoeleetrie  eonstants  (e),  and  the  elastie  eonstants  (C)  to  derive 
the  sheet  eoneentration  density  in  wurtzite  materials.  The  sheet  eoneentration  density 
cr(x)  is  given  by  ; 

\ct{x)\  =  \Pp^{Alfia,_^N)  +  Psp{Alfia,_^N)  - Psp{GaN)\ .  (3.11) 

Substituting  equation  2.9  for  the  Ppe  variable  gives: 

\cj{x)\=2  |e3;(x)-e33(x)  l  +  P^^Cxj-P^^^CO) 

V  <2q(X]  ^  [  \^C33(X]  Jj 


(3.12) 
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When  the  sheet  charge  density  <7  is  divided  by  e,  the  result  is  the  polarization 
induced  sheet  charge  density.  Figure  1 1  illustrates  the  induced  sheet  charge  as  a  function 
of  alloy  composition  (x). 
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Figure  11:  Sheet  charge  as  a  function  of  alloy  composition  [15], 


The  device  modeled  in  this  thesis  has  an  A1  alloy  content  of  0.28,  and  a  sheet  charge  of 
0.99e  cm'  was  introduced  at  the  GaN  side  of  the  heterojunction  to  simulate  the 
polarization  effects. 

Modeling  the  device  in  this  manner  resulted  in  a  more  realistic  device  model  and 
allowed  for  continued  thermal  studies.  When  modeled  with  the  interface  charge,  the 
HEMT  device  was  characterized  well,  and  with  small  discrepancies,  matched  the 
Berkeley  devices.  A  comparison  of  the  physical  device  in  relation  to  the  simulated 
device  will  be  presented  later  in  this  chapter. 
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B,  SILVACO  MODEL 

The  model  developed  in  ATLAS™  was  designed  based  on  the  work  of  Tzeng  [11], 
The  purpose  of  her  dissertation  was  to  study  noise  parameters  in  the  HEMT  deviee.  The 
Naval  Postgraduate  SehooTs  purpose  is  to  study  the  thermal  eharacteristics  of  the  device 
through  measurements  and  a  model  developed  in  Silvaco  to  determine  if  a  diamond 
substrate  will  enhance  the  thermal  characteristics  of  the  HEMT  device.  The  following 
sections  will  contain  a  discussion  of  the  modeling  steps  taken  to  simulate  the  HEMT 
device. 

1.  Model  Definition 

To  develop  the  software  model  of  the  HEMT,  ATLAS™  requires  an  input  in  the 
form  of  geometric  statements  to  define  the  dimensions  of  the  device  to  be  modeled.  The 
first  step  is  to  define  the  structure  dimensionally  and  then  assign  a  mesh  to  enclose  the 
structure.  The  purpose  of  the  mesh  is  to  define  solution  points  throughout  the  structure. 
Defining  the  mesh  correctly  is  one  of  the  most  important  factors  in  the  model 
development.  If  the  mesh  is  too  coarse,  the  solutions  will  not  be  correct,  and  if  the  mesh 
is  too  finely  defined,  the  computation  time  can  be  prohibitive.  The  rest  of  the  program 
file  is  devoted  to  modifying  statements  and  output  parameters  that  will  be  discussed 
individually  when  necessary.  The  entire  model  definition  file,  referred  to  as  a  *.m  file  is 
provided  in  Appendix  A.  Eigure  12  is  a  TONYPEOT™  structure  view  of  the  device 
modeled.  The  structure  definition  is  saved  in  ATEAS™  as  a  *.str  file.  The  graphic 
representing  the  HEMT  device  is  a  close  up  view  in  order  to  see  the  AlGaN  layer.  The 
actual  device  is  20pm  thick.  The  structure  file  view  in  TONYPEOT™  provides  a  visual 
representation  of  the  device  performance.  There  are  drop  down  menus  that  allow  for  a 
selection  of  the  desired  plot. 
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Figure  12: 2D  Representation  of  the  Modeled  AlGaN/GaN  HEMT, 

2,  Gate  Characterization 

The  characterization  of  the  gate  was  accomplished  by  measuring  the  gate  leakage 
of  the  Berkeley  device  and  then  simulating  the  device.  The  gate  on  the  HEMT  is  a 
Schottky  barrier,  and  the  characterization  was  accomplished  by  adjusting  the  gate 
workfunction  until  the  leakage  current  on  the  simulation  matched  the  measured  values. 
Figure  13  (a)  shows  the  simulated  gate  characteristic  of  the  tested  device  and  Figure  13 
(b)  shows  the  measured  characteristic  of  the  tested  device. 


TM 

Figure  13:  (a)  ATLAS  Generated  Gate  Characterization  and  (b)  Measured, 
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To  test  the  device,  the  source  and  drain  were  set  to  OV  and  the  gate  was  swept  from  -6V 
to  6V.  The  simulated  version  matches  the  tested  version  with  the  exception  of  how  low 
the  current  goes  at  Vgs=0.  The  difference  is  most  likely  due  to  the  meshing  technique 
used.  Due  to  the  magnitude  of  the  current,  is  not  a  concern  in  the  simulation.  The 
simulated  workfunction  value  of  the  Schottky  gate  that  provided  the  best  characterization 
was  4.3eV. 

3,  Heteroj unction  Bandgap  Results 

Simulating  heterojunction  devices  in  ATLAS™  is  accomplished  through 

TM  TM 

BLAZE  ,  a  compound  material  2D  simulator.  BLAZE  accounts  for  the  effects  of 
positionally  dependent  band  structure  by  modifications  to  the  charge  transport  equations 
[16].  The  difference  in  the  two  bandgaps  creates  an  abrupt  discontinuity  in  the  valence 
and  conduction  band.  ATLAS^^”,  through  BEAZE™  provides  three  methods  for  defining 
the  conduction  band  alignment  for  a  heterojunction.  The  method  utilized  in  this  model  is 

TM  TM 

the  AEIGN  parameter.  BEAZE  allows  the  difference  in  the  conduction  band  AE^,  to 

be  calculated  by  specifying  the  AEIGN  parameter  and  the  desired  value  of  offset  on  the 
material  statement.  BEAZE’^”  creates  the  desired  conduction  band  offset  by  modifying 
the  electron  affinity  of  the  material  for  which  the  parameter  is  specified.  The  AEIGN 
parameter  was  set  to  0.8  on  the  AlGaN  layer  which  allowed  80%  of  the  bandgap 
difference  to  appear  as  the  conduction  band  discontinuity. 

The  0.8  AEIGN  parameter  was  established  through  previous  HEMT  design 
attempts  [5,  7]  and  trial  and  error.  Eigure  14(a)  shows  the  valence  band  characteristic 
with  no  bias,  and  Eigure  14(b)  shows  the  conduction  band.  The  diagrams  in  Eigure  14 
are  broken  down  into  three  regions.  A  outline  was  utilized  to  extract  the  electrical 
properties  from  the  compiled  structure.  The  left  region  represents  the  gate  electrode,  the 
center  region  represents  the  AlGaN  layer,  and  the  right  region  is  the  GaN  layer.  The 
unique  bandgap  feature  of  the  HEMT  is  evident  at  the  AlGaN/GaN  interface. 
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Figure  14:  Band  diagram  of  HEMT  showing  (a)  Conduction  Band  and  (b) 

Valence  Band, 

Figure  15  is  the  full  view  of  the  bandgap  diagram  showing  the  conduction  band  in  red, 
the  Fermi  level  in  blue  and  the  valence  band  in  green.  The  bandgap  of  the  simulated 
structure  was  measured  and  found  to  be  3.42eV,  the  approximate  value  of  GaN  from 
Table  1.  The  value  of  this  figure  is  that  it  shows  the  structure  design  is  correct  and  that 
the  2DEG  is  formed  on  the  GaN  side  of  the  heterojunction. 


Bandgap  Diagram  of  NPGS  HEMT  Design 
Data  from  NPGS_SAP_g0.str 


Figure  15:  Band  Diagram  of  modeled  AlGaN/GaN  HEMT, 
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4,  Electron  Concentration  due  to  Polarization 

Figure  15  clearly  shows  the  effects  of  polarization  at  the  heterojunction  interface 


by  the  formation  of  the  triangular  shaped  quantum  well.  The  polarization  charge  causes 
the  electron  concentration  to  increase  and  a  2DEG  is  formed.  As  discussed  in  chapter  11, 
a  net  positive  polarization  charge  is  present  at  the  AlGaN/GaN  junction  in  Ga-faced 
structures.  Donor  electrons  will  attempt  to  compensate  this  charge  and  will  form  the 
2DEG  that  causes  the  conduction  band  to  drop  below  the  Eermi  level  at  the  interface. 
Eigure  16(a)  is  a  contour  view  of  the  structure  that  shows  the  electron  concentration.  The 
contour  plot  shows  that  there  is  a  distinct  2DEG  that  has  formed  just  below  the  AlGaN 
layer  (the  thin  layer  just  below  the  gate  electrode)  as  predicted.  Eigure  16(b)  shows  the 
outline  view  of  the  electron  concentration.  The  blue  line  is  the  electron  concentration  and 
the  red  line  is  the  donor  concentration.  The  increase  of  electrons  to  approximately 
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polarization  and  agrees  with  Ambacher’s  predictions[15]. 


illustrates  the  effects  of 
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Outline  Viw  of  Electtm '  ntn  ptittitmo  Donor  Concentfetion 
Data  from  NPGS_SAPj0.str 


Distance  along  line 


(a) 
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Figure  16:  Electron  Concentration  Contour  View  (a)  and  Cutline  View  (b). 
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5,  I-V  Curves 

Once  the  model  was  eomplete  and  analyzed,  a  subroutine  in  ATLAS™  was 
ereated  to  generate  the  Ids  vs.  Vds  (I-V)  eurves.  The  physieal  deviee  provided  by  Berkeley 
was  tested  utilizing  the  UTMOST™  program  and  a  Mieromanipulator  mieroscope  probe 
system.  The  data  extracted  from  the  device  was  imported  to  TONYPLOT’^”  for  analysis 
and  eomparison.  The  test  performed  on  the  model  was  performed  with  the  same  DC 
parameters  that  were  utilized  in  the  physical  device  testing.  Figure  17  is  the  I-V 
eharacteristie  of  the  HEMT  model  and  Figure  18  is  the  I-V  eharaeteristie  of  the  measured 
deviee. 


I-V  Characteristic  of  NIPGS  HEMT  Model 
Topmost  curve:  Vgs=  -1V;  step=1V 


Figure  17:  NFS  HEMT  Model  I-V  Characteristic, 
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Room  Temperature  IDS  Vs,  VDS  curve 
Topmost  curve  Vgs=-1V;  Step=1V 


Figure  18:  Measured  Berkeley  HEMT  Device  I-V  Characteristic, 

The  two  I-V  characteristic  curves  match  closely  with  a  few  exceptions.  The  topmost 
curve  on  the  modeled  device  peaks  at  about  14  mA  while  the  topmost  curve  on  the 
physical  device  peaks  at  about  12  mA.  The  topmost  curve  on  the  physical  device  has  a 
negative  slope  in  the  saturation  region  attributed  to  a  negative  resistance  due  to  the  poor 
thermal  qualities  of  sapphire.  I  was  unable  to  match  the  negative  going  slope  on  the 
modeled  device.  Although  thermal  effects  were  added  to  the  device  model,  the 
combination  used  in  this  effort  did  not  produce  the  exact  results.  The  I-V  characteristics 
match  closely  enough  to  give  me  confidence  that  the  device  was  modeled  correctly  and  it 
can  be  used  to  conduct  follow  on  thermal  measurements. 

C.  THERMAL  CHARACTERIZATION 

The  majority  of  the  modeling  efforts  in  this  thesis  have  been  dedicated  to  ensuring 
that  there  was  a  proper  HEMT  device  model  to  be  used  for  thermal  testing.  The  previous 
sections  have  demonstrated  that  the  HEMT  device  modeled  is  a  good  representation  of 
the  Berkeley  device  when  modeled  at  room  temperature.  This  section  will  analyze  the 
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model  characteristics  over  a  range  of  temperatures  and  compare  the  simulation  results  to 
the  measured  results  of  the  physical  device  over  the  same  range. 

1,  Physical  Device  Thermal  Testing 

The  thermal  testing  for  the  physical  device  was  conducted  utilizing  the 
UTMOST™  program  and  a  hot  chuck  on  the  Micromanipulator  testing  microscope.  The 
physical  device  was  tested  at  room  temperature,  100“  C,  200“  C,  and  300“  C.  Figure  19  is 
a  compilation  of  the  I-V  characteristics  over  the  range  of  temperatures. 


IDS  Reduction  as  a  function  of  Temperature 
NPGS  HEMT  Model 


Figure  19:  Topmost  Curve  (Vgs=  -1  V)  I-V  Measurement  of  Physical  Device  over 

a  range  of  temperatures. 
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Figure  20,  Figure  21,  and  Figure  22  represent  the  I-V  characteristic  of  the  room 
temperature  measurement  against  100“C,  200‘’C,  and  300“C  respectively.  The  figures 
clearly  show  the  current  degradation  due  to  the  temperature  increase. 


Figure  20:  I-V  Characteristic  Room  Temperature  vs,  100”  C, 


Figure  21:  I-V  Characteristic  Room  Temperature  vs,  200”  C, 
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Room  Temperature  vs.  300  C 
Data  from  multiple  files 


Figure  22:I-V  Characteristic  Room  Temperature  vs,  300”  C, 

Testing  the  physical  device  over  a  range  of  temperatures  served  as  the  starting  point  for 
the  thermal  model  development.  Utilizing  the  data  gained  from  the  physical  devices,  the 
simulated  device  was  modified  to  simulate  lattice  heating.  The  GIGA^m  routine  was 
utilized  in  the  thermal  testing  of  the  device  model  and  will  be  discussed  in  section  two  of 
this  chapter. 

2,  NFS  HEMT  Model  Thermal  Testing 

Introducing  the  thermal  characteristics  of  a  HEMT  device  is  accomplished 
through  the  GIGA  routine.  GIGA  extends  ATLAS  to  account  for  lattice  heat  flow 
and  general  thermal  environments  [16].  The  models  included  in  this  device  (see 
Appendix  A)  were  the  lat.temp  model  and  the  joule.heat  model.  These  two  models 
enable  ATLAS^^”  to  apply  heat  flow  equations  to  the  solution  of  the  device  parameters. 
To  simulate  the  hot  chuck  in  the  physical  device  testing,  a  thermal  contact  was  placed  in 
the  substrate  and  set  to  the  desired  temperature.  This  applied  the  same  temperature  to  the 
model  as  the  hot  chuck  applied  to  the  actual  device.  The  purpose  of  this  type  of  testing 
was  to  ensure  that  the  HEMT  model  would  react  in  the  same  manner  as  the  physical 
device  did  when  heated.  The  results,  although  not  exactly  the  same,  were  remarkably 
similar  and  indicated  that  the  modeled  device  reacted  similar  to  the  physical  device  under 
thermal  stress.  The  following  figures  are  of  the  HEMT  model  I-V  characteristics  under 
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the  same  thermal  environment  the  physical  devices  were  subjected  to.  Figure  23,  Figure 
24,  and  Figure  25  represent  the  model  I-V  characteristics  over  temperature. 


Room  Temperature  vs.  1  DO  C 
Data  from  multiple  files 


Drain  Voltage 


Figure  23: NFS  HEMT  Room  Temperature  vs,  100”  C, 


Room  Temperature  vs  200  C 
Data  from  multiple  files 


CVJ 


Figure  24:  NFS  HEMT  Room  Temperature  vs,  200”  C, 
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Room  "Temperature  vs.  300  C 
Data  from  multiple  files 


Drain  Vottage  C^J 


Figure  25:  NFS  HEMT  Room  Temperature  vs,  300”  C, 

The  I-V  characteristics  of  the  model  show  the  same  general  trend  as  the  physical  device. 
Table  4  shows  the  differences  in  the  model  and  physical  device.  The  error  analysis  is 
conducted  on  the  topmost  curve  which  represents  a  gate  voltage  of  -1  V. 


Room  Temperature 

o 

o 

o 

O 

200”  C 

300”  C 

Physical  Device 

0.012  A 

0.01  A 

0.007  A 

0.002  A 

Percent  Change 

16.6% 

41.6% 

83% 

Modeled  Device 

0.014  A 

0.009  A 

0.004  A 

0.002  A 

Percent  Change 

35.7% 

71.4% 

83% 

Table  4:  Percent  Reduction  in  Current  with  Temperature, 


Although  the  percent  change  is  approximately  double  in  the  model  design,  the  same 
general  decreasing  current  trend  is  observed  in  the  modeled  device.  The  trend  indicates 
that  the  thermal  models  incorporated  in  the  HEMT  model  are  within  a  factor  of  2  of 
degradation,  and  should  be  sufficient  to  investigate  thermal  testing.  Figure  26  illustrates 
the  decreasing  current  trend  as  the  lattice  temperature  of  the  device  is  increased  from 
room  temperature  to  300°  C. 
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Figure  26:  Decreasing  Current  trend  of  Actual  Device  and  Modeled  Device, 


An  important  measure  utilized  in  high  power  eleetronic  FET  devices  is  gate 
leakage.  When  the  device  as  previously  described,  there  is  a  significant  shift  in  the  gate 
leakage.  Figure  27  illustrates  an  approximate  increase  of  four  orders  of  magnitude  in 
gate  leakage  per  100°  C. 
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Figure  27:  NFS  HEMT  Gate  Leakage  Over  Temperature, 
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D.  FINAL  THERMAL  SIMULATION 

The  final  thermal  simulation  will  consist  of  running  the  HEMT  model  at  a  high 
DC  bias  which  will  result  in  high  temperatures  required  to  analyze  a  change  in  substrate 
materials.  To  accomplish  the  shift  in  from  sapphire  to  diamond  substrate  materials,  the 
thermal  conductivity  constant  tc.a  was  set  to  0.1  for  a  thermal  conductivity  of  10  and  0.05 
for  a  thermal  conductivity  of  20.  The  values  for  thermal  conductivity  were  taken  from 
[18],  the  SP3  corporation  website.  SP3  is  working  with  NPS  and  UC  Berkeley  to 
develop  diamond  substrates  suitable  for  AlGaN/GaN  HEMT  development.  The  device 
will  be  biased  as  described  in  the  previous  section  except  the  drain  will  be  swept  from  0 
to  55V. 

1,  Sapphire  vs.  Diamond  Substrate  Comparison 

The  point  of  testing  the  HEMT  device  on  different  substrates  was  to  determine  if 
a  change  in  substrate  material  will  reduce  the  channel  temperature.  Once  a  suitable 
HEMT  device  was  modeled,  it  was  swept  from  zero  to  55  V  on  the  drain  creating  a  high 
current  situation.  The  result  of  the  higher  drive  current  was  increased  channel/device 
temperature.  Eigure  28  shows  a  contour  plot  of  the  device  on  sapphire  and  diamond 
substrates  indicating  the  lattice  temperature  throughout  the  device.  The  thermal 
conductivity  of  diamond  ranges  between  10  and  20,  Eigure  28  represents  diamond 
thermal  conductivity  set  to  10  and  sapphire  thermal  conductivity  set  to  1.7. 
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(a)  (b) 

Figure  28:  Contour  Plot  of  (a)  Sapphire  Substrate  Material  and  (b)  Diamond 

Substrate, 
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The  plot  indicates  that  the  peak  temperature  in  the  device  changes  dramatically  with  a 
change  in  substrate  material.  Of  particular  note,  the  diamond  substrate  device  was  able  to 
dissipate  more  heat  per  unit  area,  and  the  size  of  the  “hot  spot”  was  reduced  by 
approximately  75%.  The  sapphire  substrate  device  under  the  55  V  drain  bias  reached  a 
maximum  temperature  of  157“C  while  the  diamond  substrate  device  reached  a 
temperature  of  84°  C.  Figure  30  shows  a  drain  current  of  approximately  0.06  A/mm  for 
the  diamond  substrate  device  with  a  TC  constant  of  20.  With  55  V  applied  to  the  drain  of 
the  device,  the  output  power  of  the  modeled  device  was  calculated  to  be  3.3  W/mm. 

Figure  29  is  a  plot  taken  from  a  horizontal  outline  along  the  2DEG  of  the  sapphire 
substrate  device,  the  diamond  substrate  device  with  a  thermal  conductivity  constant  of  10 
and  20.  Although  the  simulation  was  conducted  with  varying  diamond  thermal 
conductivities,  the  results  indicate  that  the  model  is  relatively  insensitive  to  the  thermal 
conductivity  change  of  diamond.  There  was  a  2%  change  in  channel  temperature  and  an 
8%  change  in  drain  current.  The  peak  temperature  of  the  device  remained  the  same  when 
the  thermal  conductivity  was  changed  from  10  to  20. 

Figure  29  better  illustrates  the  advantage  the  diamond  substrate  provides.  The 
area  under  each  of  the  thermal  curves  represents  the  average  temperature  of  the  channel. 
The  average  channel  temperature  of  the  sapphire  substrate  device  is  137°C,  the  diamond 
substrate  with  a  thermal  conductivity  parameter  of  10  reduced  the  average  channel 
temperature  to  47°C,  and  the  diamond  substrate  with  a  thermal  conductivity  parameter  of 
20  reduced  the  average  channel  temperature  to  46°C.  The  average  channel  temperature 
was  calculated  utilizing  the  integrate  function  of  TONYPLOT™.  The  simulation 
indicates  that  a  dramatic  reduction  in  lattice  temperature  and  channel  temperature  occurs 
by  changing  the  substrate  material  from  sapphire  to  diamond. 
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Channel  Temperature  Thermal  Plot 
NPGS  HEMT  Model 


Figure  29:  Thermal  Horizontal  Cutline  of  2DEG  with  Diamond  Thermal 
Conductivity  set  to  10  and  20,  Sapphire  (hrown)  Diamond  10  (green)  and 

Diamond  20  (red). 


Figure  30  shows  the  device  drain  current  utilizing  the  three  different  substrates. 


Drain  Current  Plot 
NPGS  HEMT  Model 


Figure  30:  Drain  Current  as  a  Function  of  Substrate  Material, 
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The  final  parameter  to  be  analyzed  is  the  change  in  the  gate  leakage  characteristic  with  a 
change  in  substrate  material.  Figure  31  illustrates  the  gate  leakage  characteristic  as  a 
function  of  substrate  material. 
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Figure  31:  Gate  Leakage  as  a  Function  of  Substrate  Material, 


There  is  a  reduction  of  five  orders  of  magnitude  in  the  amount  of  gate  leakage  on  the 
diamond  substrate  device  compared  to  the  sapphire  substrate  device.  Every  measure  of 
effectiveness  of  a  HEMT  device  improved  with  a  change  in  substrate  material.  The 
simulation  indicates  that  HEMT  devices  grown  on  diamond  substrates  will  improve  the 
overall  effectiveness  of  the  device  when  high  power  applications  are  required.  Although 
it  is  beyond  the  scope  of  this  thesis,  it  is  logical  to  assume  that  the  reliability  of  the 
devices  will  be  enhanced  if  the  temperature  can  be  reduced. 

Table  5  shows  the  results  of  the  various  parameters  tested.  The  change  was 
calculated  with  the  diamond  TC=10  parameter. 
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Diamond  TC=10 

Diamond  TC=20 

Sapphire 

Change 

Peak  Temperature 

84  C 

84  C 

150  C 

56% 

Channel  Temperature 

47  C 

46  C 

137  C 

34% 

Gate  Leakage 

7.00E-20 

2.00E-20 

1.00E-15 

1E5 

Drain  Current 

.57  A 

.62  A 

.55  A 

0.04% 

Table  5:  Summary  of  Testing  Results, 


E,  CHAPTER  CONCLUSION 

This  chapter  introduced  the  modeling  techniques  utilized  to  simulate  the 
AlGaN/GaN  HEMT  deviee.  To  ensure  proper  thermal  modeling  a  HEMT  model  was 
developed  and  tested  both  at  room  temperature  and  at  various  temperatures.  The  results 
of  the  model  were  compared  with  the  physieal  device  for  aecuraey.  Once  a  suitable 
device  strueture  was  developed,  it  was  subjeeted  to  a  high  bias  condition  to  generate 
internal  heating.  The  device  parameters  were  then  extracted  from  a  sapphire  based 
substrate  and  then  a  diamond  substrate  to  determine  if  the  substrate  material  enhanced  the 
thermal  capacity  of  the  device.  It  was  shown  that  a  diamond  substrate  has  a  dramatie 
effect  on  channel  temperature.  There  was  a  54%  reduction  in  peak  temperature  and  a 
34%  reduetion  in  the  channel  temperature.  The  following  chapter  briefly  summarizes  the 
work  eompleted  on  this  thesis  and  will  make  recommendations  for  follow-on  simulations. 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  AlGaN/GaN  HEMT  model  developed  in  this  thesis  provides  a  representation 
of  an  actual  HEMT  device  that  was  tested  over  a  range  of  temperatures  from  21^C  to 
300°C  and  validated  against  measured  data.  The  thermal  effects  that  occur  by  device  self 
heating  effects  were  correctly  modeled  through  the  Silvaco’s  GIGA™  routine.  The 
HEMT  model  was  able  to  duplicate  the  actual  transistor  operation  both  electrically  and 
thermally.  The  I-V  curves  generated  by  the  software  model  represented  the  actual  device 
I-V  curves  well  with  only  a  13%  difference  at  the  topmost  I-V  curve. 

The  difficult  part  of  this  thesis  was  ensuring  that  a  proper  HEMT  device  was 
modeled.  A  tremendous  effort  was  made  ensuring  all  operational  aspects  of  the  HEMT 
device  were  addressed.  The  focus  of  the  thesis  however,  was  to  examine  the  thermal 
effects  of  the  device  when  the  substrate  material  was  changed  from  sapphire  to  diamond. 
Changing  the  substrate  and  examining  the  thermal  properties  was  accomplished  by 
changing  the  material  name  in  the  ATLAS™  file  and  adjusting  the  thermal  conductivity 
constants.  Sapphire’s  thermal  conductivity  was  set  to  1.7  W/cm-K,  and  diamond’s 
thermal  conductivity  was  set  to  both  10  and  20  W/cm-K  to  test  over  the  range  given  by 
SP3  incorporated.  The  effects  of  the  substrate  change  were  drastic.  There  was  a  1/2 
reduction  in  the  peak  temperature  of  the  device  when  diamond  was  used  see  Table  5. 
The  average  temperature  in  the  channel  was  reduced  by  1/3.  Thermal  analysis  is 
important  in  the  calculation  of  device  lifetime.  One  approach  to  predicting  device 
median  lifetimes  is  the  Arrhenius  model  for  temperature  dependence  [4]  and  is  given  by: 

tm('r)  =  Ae^^  (4.1) 

where  EA  is  the  activation  energy  (eV),  A  is  a  constant  acceleration  factor 
(dimensionless)  and  k  is  the  Boltzmann’s  constant  (eV/T).  Activation  energies  range 
from  0.3  eV  to  1.7  eV  for  GaN  HEMTs  [7].  Utilizing  an  acceleration  factor  of  1  and 
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activation  energy  of  0.3  eV  and  1.7  eV,  the  median  lifetime  of  a  HEMT  deviee  can 
theoretically  be  increased  to  a  range  of  5  to  le6  times  the  median  lifetime  of  a  sapphire 
substrate  HEMT. 

Initially,  1  attempted  to  utilize  the  new  polarization  models  available  in 
ATEAS'^’’^.  Many  attempts  to  model  the  polarization  effects  by  utilizing  the  PSP  and 
POEARIZATION  statements  failed.  1  was  unable  to  get  a  suitable  bangap  diagram  or 
2DEG  utilizing  the  new  BLAZE™  models.  The  interface  statement  on  the  AlGaN  side 
of  the  heterojunetion  provided  the  best  method  to  simulate  the  polarization  effeets 
experienced  in  wurtzite  materials.  The  interface  command  places  a  sheet  eharge  at  the 
interface  and  provides  the  necessary  eleetron  concentration  to  form  the  2DEG. 

The  model  designed  in  this  thesis  proved  to  be  relatively  insensitive  to  the  thermal 
conductivity  constants  given  for  diamond.  There  was  very  little  change  in  any  of  the 
parameters  analyzed  when  the  constants  were  varied  through  the  published  values. 
Although  many  factors  in  the  ATLAS’’’’^  code  had  to  be  optimized  to  get  a  representative 
HEMT  deviee,  one  of  the  enabling  parameters  to  simulate  a  geometrically  correct  device 
was  the  Albreeht  model.  Albreeht  et  al.  [19]  developed  a  set  of  transport  parameters  for 
devices  with  eleetron  conduction  GaN  (position  of  the  2DEG).  The  Albreeht  parameter 
models  the  carrier  mobility  as  a  function  of  doping  and  lattice  temperature  which  are 
critical  parameters  in  the  deviee  modeled  for  this  thesis.  In  [19]  Albreeht  provides  initial 
experimental  values  for  the  variables  available  in  the  ATLAS™  model.  These  values 
were  implemented  in  the  HEMT  model  and  then  optimized  to  match  the  measured  device 
charaeteristies.  The  variables  in  the  Albreeht  model  proved  to  be  the  most  sensitive 
parameters  affeeting  the  I-V  curves  and  further  analysis  and  optimization  efforts  should 
provide  a  better  I-V  representation. 


B.  RECOMMENDATIONS 

Eollow-on  efforts  should  examine  the  possibility  of  implementing  the  polarization 
and  piezoelectric  models  available  in  ATLAS’’’’^  to  determine  if  a  suitable  model  can  be 
developed  utilizing  the  new  polarization  models. 
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The  parameters  examined  in  this  thesis  were  DC  measurements  only.  One  of  the 
advantages  of  a  HEMT  is  that  it  has  the  ability  to  operate  at  high  frequeneies.  It  would 
be  beneficial  to  examine  the  RF  characteristics  of  the  model  developed  in  this  thesis  to 
see  if  the  modeled  device  can  accurately  simulate  the  RF  parameters  of  an  actual  HFMT 
device.  Once  the  basic  characteristics  are  extracted,  the  device  should  be  simulated  at 
high-frequency,  high-power  conditions  to  determine  if  the  benefits  of  a  diamond  substrate 
are  as  dramatic  in  the  AC  domain. 

The  HFMT  design  presented  utilized  an  18  pm  thick  sapphire  or  diamond 
substrate.  The  purpose  of  simulating  the  model  with  a  thick  substrate  was  to  analyze  the 
thermal  parameters  when  the  device  was  not  saturated  thermally.  The  assumption  in  this 
model  is  that  the  device  is  capable  of  dissipating  all  of  the  heat  applied.  In  reality,  the 
substrate  thickness  will  be  much  smaller  and  a  more  comprehensive  heat  flow  study 
should  be  conducted  in  order  to  simulate  a  more  realistic  environment. 

The  two  dimensional  analysis  of  the  device  provided  and  adequate  simulation  of 
the  HFMT  device,  however  Silvaco  has  a  three  dimensional  capability  that  can  provide  a 
more  comprehensive  analysis  of  the  device  operational  parameters.  The  model 
developed  in  this  thesis  can  be  expanded  to  a  three  dimensional  model  in  order  to  re¬ 
examine  the  electrical  and  thermal  properties  of  the  device.  An  initial  attempt  was  made 
to  extend  the  model  to  three  dimensions,  but  I  was  unable  to  get  the  model  to  converge  at 
a  high  bias  condition.  The  probable  cause  of  the  failure  was  meshing  techniques.  Great 
care  must  be  taken  to  properly  mesh  the  device  for  simulation.  In  three  dimensions,  the 
number  of  nodes  required  for  a  solution  increases  dramatically.  If  there  are  too  many 
nodes,  the  computation  time  is  prohibitively  long,  and  if  the  mesh  is  reduced  too  much, 
the  results  will  not  be  accurate. 

The  physical  device  grown  on  a  sapphire  substrate  tested  for  this  thesis  was 
designed  and  fabricated  by  S.  Tzeng  at  UC  Berkeley.  The  Weber  material  science  group 
at  UC  Berkeley,  NFS,  and  SP3  are  currently  working  to  grow  AlGaN/GaN  HEMT 
devices  on  diamond  substrates.  Once  a  viable  diamond  substrate  device  is  created,  a 
comprehensive  test  should  be  preformed  on  the  device  to  see  if  the  ATFAS™  model 
accurately  predicts  the  thermal  effects  of  the  diamond  substrate  device.  Once  the 
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diamond  HEMTs  are  available  and  tested  to  ensure  model  aceuraey,  eleetrical  defeet 
densities  could  be  modeled  to  increase  with  temperature.  A  study  could  be  performed  to 
predict  degradation  due  to  temperature  and  time.  The  initial  results  gained  from  this 
thesis  indicate  that  AlGaN/GaN  HEMTs  grown  on  a  diamond  substrate  have  the 
possibility  to  be  the  enabling  technology  for  high-power  high-frequency  operations. 
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APPENDIX:  ATLAS^m  INPUT  DECK 


###################################################################### 

###################################################################### 

#LT  Scott  Newham 

#NPS  AlGaN/GaN  HEMT  Design 

# 

#267  Angstrom  AlGaN  layer 

#  Schottky  WF=4 . 3 

#  lel5  doping  GaN  AlGaN=lelO 
#Interface  charge=. 99el3 

#  no  doping  in  AlGaN 

#  Albrct  mobility  model  an . albrct=5e-3  bn . albrct=3e-3  cn . albrct=2e-3 
####################################################################### 
go  atlas 

set  devthk=20 

set  sourcegatespace=2 

set  gatedrainspace=3 

set  devwidth=4+$sourcegatespace+$gatedrainspace 

set  reg3xmin=$devwidth-l 

set  reg4xmin=l+$sourcegatespace 

set  reg4xmax=2+$reg4xmin 

set  reg5xmax=$devwidth- . 5 

set  drainxmin=$devwidth-l 

set  drainxmax=$devwidth- . 6 

set  ifchargemax=$devwidth-l . 05 

set  WF=4.3 

#  Sweep  increment 
set  vstart  =  0 
set  vstop  =  7 

set  vine  =  .25 

####################  Mesh  Construction  ############# 

####Mesing  and  width  =  100  microns################# 

mesh  auto  width=100 

x.m  1=0 . 0  s=0 . 1 

x.m  1=0 . 5  s=0 . 1 

x.m  1=1 . 0  s=0 . 1 

x.m  l=$reg4xmin  s=0 . 1 

x.m  l=$reg4xmax  s=0 . 1 

x.m  l=$reg3xmin  s=0 . 1 

x.m  l=$reg5xmax  s=0 . 1 

x. m  l=$devwidth  s=0 . 1 

y. m  1=0  s=0.5 

y .m  1=0 . 4  s=0 . 1 
y.m  1=0.49  s=0.005 
y.m  1=0.5  s=0.005 
y.m  1=0.525  s=0.005 
y.m  1=0.5267  s=0. 00005 
y.m  1=0.5269  s=0. 00005 
y.m  1=0.5270  s=0.005 
y.m  1=1 . 0  s=0 . 1 
y.m  1=2  s=0 . 5 
y.m  l=$devthk  s=2 . 0 

########################  Region  Definitions  ################ 

region  num=l  mat=AlGaN  donors=lel0  x.comp=0.28  x.min=0.5  x . max=$reg5xmax 
y.min=0.5  y.max=0.5267 

region  num=2  mat=oxide  x.min=0  x . max=$devwidth  y.min=0  y.max=0.5 
region  num=3  mat=GaN  donors=lel5  x.min=0.5  x . max=$reg5xmax  y.min=0.5267 
y.max=0 .5269 
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region  nuin=4  inat=GaN  donors=lel5  x.min=0.5  x . max=$reg5xmax  y . inin=0 . 52 69 
y .max=l . 0 

region  nuin=5  inat=GaN  donors=lel5  x.min=0  x . max=$devwidth  y.min=l  y.max=2 
region  num=6  mat=Sapphire  x.min=0  x .max=$devwidth  y.min=2  y .max=$devthk; 
elec  num=l  name=source  x.min=0  x.max=0.5  y.min=0.5  y.max=l 
elec  num=2  naine=drain  x .min=$drainxmax  x . inax=$devwidth  y.inin=0.5  y.max=l 
elec  num=3  name=gate  x.min=$reg4xinin  x . max=$reg4xmax  y.min=0  y.max=0.5 
elec  num=4  substrate 

##############  Modifying  Statements  ################## 

#######Inserts  polarization  effects####### 
interface  charge=0 . 99el3  y.min=.52  y.max=.5267  s.s 

material  mat=AlGaN  align=.8 

mobility  albrct.n  an . albrct=5e-3  bn . albrct=3e-3  cn . albrct=2e-3  vsatn=2e7 

material  mat=GaN  tcon.polyn 

material  mat=AlGaN  tcon.polyn 

material  mat=Sapphire  tcon.polyn 

models  k.p  print  lat.temp  joule. heat  srh  albrct 

contact  name=gate  work;=$WF 

#thermcontact  num=l  x.min=0  x .max=$devwidth  y.min=2  y.max=5  temp=300  ^boundary 
alpha=l . 7 

thermcontact  num=l  x.min=0  x . max=$devwidth  y.min=14  y .max=$devthk  temp=300 

^boundary  alpha=l . 7 

output  con. band  val.band  charge 

method  gumits=300  clim.dd=le5  autonr  block  carr=l 
######################  Output  Statements  ########## 

#  idvd  curves 
solve 

save  outf =NPGS_SAP_gO . str 
log  outf=NPGS_SAP_d_0 . log 

solve  name=drain  vdrain=$vstart  vf inal=$vstop  vstep=$vinc 
save  outf =NPGS_2D_SAP_g0d7 . str 
log  off 

tonyplot  NPGS_SAP_d_0.log  -set  IDVD. set 
solve  vdrain=0  vgate=  -5 
save  outf=solve_vgate-5 . str 
log  outf=NPGS_SAP_d_l . log 

solve  name=drain  vdrain=$vstart  vf inal=$vstop  vstep=$vinc 
log  off 

solve  vdrain=0  vgate=  -4 
save  outf=solve_vgate-4 . str 
log  outf=NPGS_SAP_d_2 . log 

solve  name=drain  vdrain=$vstart  vf inal=$vstop  vstep=$vinc 
log  off 

solve  vdrain=0  vgate=  -3 
save  outf=solve_vgate-3 . str 
log  outf=NPGS_SAP_d_3 . log 

solve  name=drain  vdrain=$vstart  vf inal=$vstop  vstep=$vinc 
log  off 

solve  vdrain=0  vgate=  -2 
save  outf=solve_vgate-2 . str 
log  outf=NPGS_SAP_d_4 . log 

solve  name=drain  vdrain=$vstart  vf inal=$vstop  vstep=$vinc 
log  off 

solve  vdrain=0  vgate=  -1 
save  outf=solve_vgate-l . str 
log  outf=NPGS_SAP_d_5 . log 

solve  name=drain  vdrain=$vstart  vf inal=$vstop  vstep=$vinc 
log  off 

#solve  vdrain=0  vgate=  0 
#save  outf=solve_vgate_0 . str 
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#log  outf=NPGS_SAP_d_6 . log 

#solve  naine=drain  vdrain=$vstart  vf inal=$vstop  vstep=$vinc 
#log  off 

#solve  vdrain=0  vgate=  1 
#save  outf=solve_vgate_l . str 
#log  outf =NPGS_SAP_d_7 . log 

#solve  naine=drain  vdrain=$vstart  vf inal=$vstop  vstep=$vinc 
#log  off 

tonyplot  NPGS_SAP_d_0.log  -overlay  NPGS_SAP_d_l.log  -overlay  NPGS_SAP_d_2.log  - 
overlay  NPGS_SAP_d_3.log  -overlay  NPGS_SAP_d_4.log  -overlay  NPGS_SAP_d_5.log  -\ 
overlay  NPGS_SAP_d_6.log  -overlay  NPGS_SAP_d_7.log  -set  IDVD.set 
#  idvg  curve 

#log  outf =NPGS_GaNsapp-gate_sweep . log 
#solve  vdrain=0 

#solve  naine=gate  vgate=0  vfinal=-4.2  vstep=.6 
quit 
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